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ABSTRACT 

Based on numerical simulations of gas-rich galaxy mergers, we discuss a model in which quasar activity is 
tied to the self-regulated growth of supermassive black holes in galaxies. The nuclear inflow of gas attending a 
galaxy collision triggers a starburst and feeds black hole growth, but for most of the duration of the starburst, the 
black hole is "buried", being heavily obscured by surrounding gas and dust, limiting the visibility of the quasar, 
especially at optical and ultraviolet wavelengths. As the black hole grows, feedback energy from accretion heats 
the gas and eventually expels it in a powerful wind, leaving behind a "dead quasar". In between the buried and 
dead phases, there is a window in time during which the galaxy would be seen as a luminous quasar. Because 
the black hole mass, radiative output, and distribution of obscuring gas and dust all evolve strongly with time, 
the duration of this phase of observable quasar activity depends on both the waveband and imposed luminosity 
threshold. We determine the observed and intrinsic lifetimes as a function of luminosity and frequency, and 
calculate observable lifetimes ~ lOMyr for bright quasars in the optical B-band, in good agreement with 
empirical estimates and much smaller than our estimated black hole growth timescales ~ lOOMyr, naturally 
producing a substantial population of buried quasars. However, the observed and intrinsic energy outputs 
converge in the IR and hard X-ray bands as attenuation becomes weaker and chances of observation greatly 
increase. We also obtain the distribution of column densities along sightlines in which the quasar is seen above a 
given luminosity, and find that our result agrees remarkably well with observed estimates of the column density 
distribution from the SDSS for the appropriate luminosity thresholds. Our model reproduces a wide range of 
quasar phenomena, including observed quasar lifetimes, intrinsic lifetimes, column density distributions, and 
differences between optical and X-ray samples, having properties consistent with observations across more 
than five orders of magnitude in bolometric luminosity from 10 9 - 10 14 L Q (-17 < Mb < -30). 

Subject headings: quasars: general — galaxies: nuclei — galaxies: active — galaxies: evolution — cosmology: 
theory 



1. INTRODUCTION 

More than 40 years have passed since the recognition that 
quasars are at cosmological distances (Schmidt 1963; Green- 
stein & Mathews 1963; Mathews & Sandage 1963) and hence 
must be powerful energy sources. A combination of argu- 
ments based on time variability and energetics strongly sup- 
ports the view that the activity is produced by the accretion of 
gas onto supermassive black holes in the centers of galaxies 
(e.g. Salpeter 1964; Zel'dovich & Novikov 1964; Lynden- 
Bell 1969). However, the mechanism that provides the trigger 
to fuel quasars remains uncertain. Recent discoveries of cor- 
relations between masses of black holes in nearby galaxies 
and either the mass (Magorrian et al. 1998) or velocity dis- 
persion (i.e. the Mbh-c relation: Ferrarese & Merritt 2000; 
Gebhardt et al. 2000) of spheroids demonstrate a fundamen- 
tal link between the growth of supermassive black holes and 
galaxy formation. However, theoretical evidence connecting 
the origin of supermassive black holes and quasars to galaxy 
evolution remains elusive, owing to the complexity of the un- 
derlying physics and dynamics. In particular, up to now there 
has been no comprehensive model to explain the origin and 
fueling of quasars, or their lifetimes, obscuration, demograph- 
ics, self-regulation and termination, and dependence on host 
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galaxy properties. 

Observations of the nearby Universe suggest that rapid 
black hole growth may be related to massive flows of gas into 
the centers of galaxies. Infrared (IR) luminous galaxies repre- 
sent the dominant population of objects above 10 n Lo locally 
and it is believed that much of their IR emission is powered by 
dust reprocessing of radiation from intense nuclear starbursts 
(e.g. Soifer et al. 1984a,b; Sanders et al. 1986, 1988a,b; 
for a review, see e.g. Soifer et al. 1987). At the highest lu- 
minosities above 10 12 L Q , characteristic of ultraluminous in- 
frared galaxies (ULIRGs), nearly all the galaxies appear to 
be in advanced stages of merging (e.g. Allen et al. 1985; 
Joseph & Wright 1985; Armus et al. 1987; Kleinmann et 
al. 1988; for reviews, see Sanders & Mirabel 1996 and Jo- 
gee 2004) and CO observations show that they contain large 
quantities of gas in their nuclei (e.g. Scoville et al. 1986; Sar- 
gent et al. 1987, 1989). Some ULIRGs exhibit "warm" IR 
spectral energy distributions, perhaps indicative of a buried 
quasar (e.g. Sanders et al. 1988c). This fact, together with 
the overlap between bolometric luminosities of ULIRGs and 
quasars, led Sanders et al. (1988a) to propose that quasars are 
the descendents of an infrared luminous phase of galaxy evo- 
lution caused by mergers. This scenario is supported by re- 
cent X-ray observations which have revealed the presence of 
two non-thermal poi nt sources n ear the center of the ULIRG 
NGC6240 ( Komossa et a 1. 2003), which are most naturally in- 
terpreted as accreting supermassive black holes that are heav- 
ily obscured at visual wavelengths (e.g. Gerssen et al. 2004; 
Max et al. 2005). 

Hydrodynamic simulations have shown that gas inflows are 
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produced by strong gravitational torques on the gas through 
tidal forces duri ng mergers involving gas-rich gala xies (e.g., 
Hernquist 19891 iBarnes & Hernom^ll99"ll fT996h . The re- 
sulting dense concentrations of gas in the inner regions of 
the remnant can sustain star formation at rates high enough 
and for sufficiently long timescales to account for many 
properties of ULIRGs (e. g., iMihos & Hernauisllfl99H fl996l 
iHernauist & MihoslfT995h . However, these authors were not 
able to directly explore the relationship between the gas in- 
flows and quasar activity because their simulations did not in- 
clude a model for black hole growth and the impact of strong 
feedback from either star formation or quasars. 

One of the most fundamental parameters of black hole 
growth is the quasar lifetime, fg, which sets the timescale 
for the most luminous phase of the activity. Observa- 
tions generally constrain quasar lifetimes to the range ?g w 
10 6 - 10 8 yr (for a review, see Martini 2004J). These 
estimates are primarily based on demographic or inte- 
gral arguments which combine observations of the present- 
day population of supermassive black holes and accretion 
by the high-redshift quasar popula t ion ("e.g.. [Soltan l [19821: 
t Nataraian. & Reesl Tl998t lYu & Tremamel hooi 



Yu & Lu 2004; Haiman et al. 2004), or incorporate quasars 
into models of galaxy evolution ("e.g..|Ka uffmann & Haehnelt 
[ 20001 IWvithe & Loebl 120021 iDi Matteo et all 120031 120041 
Grana to et alJ l2004) or reionization of Hell (e.g. Sokasian 
et al. 2002, 2003). Results from clustering in quasar surveys 
(e.gJPorciani. Magliocchetti. & Norberg|2004tlGrazian et alJ 



2004), the proximity effect in the Lya forest (Bajtlik, Dun- 
can & Ostriker 1988; Haiman & Cen 2002; Yu & Lu 
2005; although such results may still be inconclusive, e.g. 
Croft 2004), and the transverse proximity effect in He II 
( Jako bsen et all2003l) similarly suggest lifetimes tg ~ 10 7 yr. 

It is not immediately obvious how quasar lifetimes tQ ~ 
10 7 yr can be explained in the context of gas inflows triggered 
by galaxy mergers. The full duration of the starburst phase is 
set by the timescale during which strong gas inflows are ex- 
cited, which in turn is determined by the time when significant 
gravitational t orques are exerted on the gas. T he numerical 
simulations of Mihos & Hernquist ( 1994, 1996) showed that 
this occurs for ^ 2 x 10 8 yr, much shorter than typical m erger 
times cales of > 10 9 yr JBarnesll 19881 1 19921 iHernquistll 19921 
119931) . but in good agreement with observed estimates for the 
gas depletion time in the central regions of ULIRGs (for a 
discussion, see e.g. Barnes & Hernquist 1992). If gas-rich 
mergers are indeed responsible for the origin of quasar activ- 
ity, as suggested through studies of ULIRGs (e.g. Sanders 
et al. 1988a) and more directly from observations of quasar 
hosts (e.g. Stockton 1978; Heckman et al. 1984; Stockton 
& MacKenty 1987; Stockton & Ridgway 1991; Hutchings & 
Neff 1992; Bahcall et al. 1994, 1995; Canalizo & Stockton 
2001), the factor ^10 difference between the gas consump- 
tion time in ULIRGs and quasar lifetimes must be reconciled. 

Semi-analytical models of supermassive black hole evo- 
lution and its correlation with galaxy structure suggest that, 
beyond a certain threshold, feedback en ergy expels nearby 
gas and shuts down the accr etion p hase (|Silk & Reesl 1 19981 
Fabian 1999; Wyithe &Tol^l2003l However, these calcu- 
lations provide only limits, and in neglecting the dynamics 
of quasar evolution they do not predict time-dependent ef- 
fects such as the characteristic lifetime of the accretion phase 
prior to its self-termination, the fueling rates for black hole 
accretion, the obscuration of central sources, or the quasar 



light curve. These quantities are instead usually taken to be 
independent input parameters of the model. For example, 
the lifetime is either adopted from observational estimates or 
assumed to be similar to a characteristic timescale such as 
the dynamical time of the host galaxy disk or the e-folding 
time for Eddington-limited black hole growth t$ = M^n/M = 
4.5 x 10 7 /(e r /0.1)yr for accretion with radiative efficiency 
e,. = L/Mc 2 ~ 0.1 and I = L/L Edd < 1 (Salpe terl 19641 

Efforts to model quasar accr etion and feedback in a more 
self-consistent manner (e.g., iCiotti &OstrikeH[T997i l200~fl) 
have considered the time-dependent nature of gas inflows 
feeding accretion and the impact of radiative heating from 
quasar feedback in a hydrodynamical context. Although gen- 
erally restricted to one-dimensional "toy models" which do 
not include galaxy-galaxy interactions, such modeling has 
made important progress in predicting the char acteristic duty 
cycles of quasars, in good agreement with e.g., Haiman et al. 
( 2004), and in demonstrating the importance of radiative feed- 
back on the surrounding medium and quasar evolutio n (for a 
review , see lOstriker & Ciottil f2005). Further. fSazonov et alJ 
( 120051) showed that self regulation by radiative feedback from 
quasars in this modeling can expel a significant quantity of 
gas and leave a remnant with properties similar to observed 
ellipticals and black hole mass corre sponding to the observed 
Mbh-c relation. Similarly Kawata & Gibson (2005) showed 
that incorporating a simple prescription for radiative heating 
by quasars in the context of cosmological simulations yields 
populations of elliptical galaxies consistent with observed op- 
tical and X-ray properties. 

Recently, Springel et al. (2005b) have developed a method- 
ology for incorporating black hole growth and feedback into 
hydrodynamical simulations of galaxy mergers that includes a 
multiphase model for star formation and pressurizatio n of the 
inters tellar gas by supernova feedback ( Springel & Hernauist 
120031) . Using this approach, we have begun to explore the im- 
pact of these processes on galaxy formation and evolution. Di 
Matteo et al. (2005) and Springel et al. (2005b) have shown 
that the gas inflows produced by gravitational torques during 
a merger both trigger starbursts (as in the earlier simulations 
of e.g. Mihos & Hernquist 1994, 1996) and fuel rapid black 
hole growth. The growth of the black hole is determined by 
the gas supply and terminates abruptly when significant gas 
is expelled owing to the coupling between feedback energy 
from black hole accretion and the surrounding gas. Eventu- 
ally, as the gas is heated and driven out, the remnant is no 
longer active because the black hole does not accrete at a high 
rate, leaving a dead quasar in an ordinary galaxy. The self- 
regulated nature of the black hole growth explains observed 
correlations between black hole ma ss and properties of nor- 
mal galaxies (Di M atteo et alj|2005| ). as well as the color dis- 
tribution of ellipticals ( Spring el et al. Il2005al) . These results 
lend support to the view that mergers have played an impor- 
tant role in structuring galaxies, as advocated by e.g. Toomre 
& Toomre (1 972) and Toomre (1977 ). 

Moreover, Di Matt eo et al.l (|2005) showed that the dynam- 
ics of the inflowing gas and its response to the self-regulated 
growth of the black hole yield a timescale for the strong ac- 
cretion phase ~2x 10 s yr, comparable to the full duration 
of the starburst. During much of this period, the bolometric 
luminosity of the black hole would exceed the threshold to 
be classified as a quasar, implying a timescale ~ 10 8 yr (for a 
Milky Way mass system) for the intrinsic quasar phase of the 
black hole. 
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In lHopkins et al.l ( 12005 al) . we employed models for obscu- 
ration by gas and dust to show that for most of the accretion 
lifetime the quasars would be buried; i.e. heavily obscured 
by the large gas density powering accretion. Eventually, feed- 
back from the accretion energy drives away the gas, creating a 
brief window in which the central object would be observable 
as an optical quasar, until accretion levels drop below quasar 
thresholds. By calculating the effects of obscuration from 
the simulation of a merger of gas-rich galaxies, we showed 
that this determines an observable lifetime ~ 10 7 yr, in good 
agreement with observations. 

Here, we extend and further devel op implications of the 
model for quasar evolution proposed in Hop kins et al.l f2005a) 
by analyzing a series of hydrodynamical simulations of 
galaxy mergers where we vary the masses of the progenitor 
galaxies, so that they have virial velocities between 80 and 
320km s" 1 . In §2 we describe our series of merger simula- 
tions and the calculation of column densities and obscuration, 
and compare the results obtained with different column den- 
sity calculations. In §3 we examine the frequency dependence 
of our results and quantify the differences in lifetimes across 
the optical and X-ray bands. In §4 we examine the typical col- 
umn densities along lines-of-sight to the simulated quasars, 
and compare these to observed distributions from both opti- 
cal and X-ray surveys. In §5 we show the simulation results 
for our merger series and compare the observed and intrinsic 
lifetimes as a function of luminosity for all the simulations. 
Finally, in §6 we discuss our results and the implications of 
our model for a broad range of quasar studies. 

2. THE MODEL: CALCULATIONS 
2.1. The Simulations 

The simulations were performed using the GADGET- 
2 code, a new version of the parallel T reeSPH code 
GADGET (iSpringel. Yoshida. & White! 120011). It uses an 
entropy-conserving formulation ( Springel & Hernauist 2002) 
of smoothed particle hydrodynamics (SPH), and includes ra- 
diative cooling, heating by a UV background (as described in 
e.g. Katz et al. 1996; Dave et al. 1999), and a sub-resolution 
model of a multiphase interstellar medium (ISM) to describe 
star fo rmation and supernova feedback ( Spri ngel & Hernauist! 
2003). This sub-resolution model provides an effective equa- 
tion of state for star-forming gas which includes pressure 
feedback from supernova heating, and allows us to stably 
evolve even massive pure gas disks (see, e.g. Springel et al. 
2005b; Robertson et al. 2004). The methodology of accre- 
tion, feedback , and galaxy generation is described in detail in 
Spri ngel et aD ( l2005bl) . 

In our approach, supermassive black holes (BHs) are repre- 
sented by "sink" particles that accrete gas from their local en- 
vironment, with an accretion rate M estimated from the local 
gas density and sound speed using a Bondi-Hoyle-Lyttleton 
parameterization with an imposed upper limit equal to the Ed- 
dington rate. The bolometric luminosity of the BH particle is 
then Lboi = £ r Mc 2 , where e r = 0.1 is the radiative efficiency. 
We further assume that a small fraction (5%) of Lboi couples 
dynamically to the surrounding gas, and that this feedback 
is injected into the gas as thermal ene rgy. This fract ion is a 
free parameter, determined in Di Matte o et alJ (|2005) by fit- 
ting to the Mbh-c relation. We do not attempt to resolve 
the small-scale accretion dynamics near the black hole, but 
instead assume that the time-averaged accretion can be esti- 
mated from the gas properties on the scale of our spatial reso- 



TABLE 1 
Simulation Parameters 



Simulation 


Vvirtkms- 1 ) 


M vir (M ) 


M B H. final (M ) 


Al 


80 


1.7 X 10" 


7 x 10 6 


A2 


113 


4.8 X 10" 


3 x 10 7 


A3 


160 


1.4 X 10 12 


3 x 10 8 


A4 


226 


3.8 x 10 12 


7 x 10 8 


A5 


320 


1.1 x 10 13 


2x 10 9 



lution (< 50 pc). 

In what follows, we analyze five simulations of colliding 
disk galaxies, which form a family of structurally similar 
models with different virial velocity and mass. In each sim- 
ulation, we generate two stable, isolated d isk galaxies, each 
with an extended dark matter halo with a Hern auisl i 19901) 
profile, motivated by cosmological simulations (e.g. Navarro 
et al. 1996; Busha et al. 2004) and observations of halo 
properties (e.g. Rines et al. 2002, 2002, 2003, 2004), an 
exponential gas disk, and a b ulge. The simu l ations follow 
the series described in detail in Springel et al. (2005a), with 
the parameters listed in Tabled We denote the simulations 
Al, A2, A3, A4, and A5, with increasing virial velocities 
of V vir = 80, 113, 160,226,and 320^8"', respectively. Note 
that the self-similarity of these models is broken by the scale- 
dependent phys ics of cooling, star f ormation, and black hole 
accretion. In Hopki nTet alJ ll2005al) . we describe our analy- 
sis of simulation A3, a fiducial choice with a rotation curve 
and mass similar to the Milky Way. The galaxies have mass 
M v j r = V v 3 ir /(10G//oX with the baryonic disk having mass frac- 
tion f7?d = 0.041, the bulge ntb = 0.0136, and the rest of the 
mass in dark matter with a concentration parameter 9.0. The 
disk-scale length is computed based on an assumed spin pa- 
rameter A = 0.033, and the scale-length of the bulge is set to 
0.2 times this. We begin our simulation with pure gas disks, 
which may better correspond to the high-redshift galaxies in 
which most quasars are observed. 

Each galaxy is initially composed of 168000 dark matter 
halo particles, 8000 bulge particles, 24000 gaseous disk par- 
ticles, and one BH particle, with a small initial seed mass of 
1O 5 M . Given these choices, the dark matter, gas, and star 
particles are all of roughly equal mass, and central cusps in the 
dark matter and bulge profiles are reasonably well resolved 
(see Fig 2. in Springel et al. 2005b). The galaxies are then set 
to collide in a prograde encounter with zero orbital energy and 
a small pericenter separation (7.1 kpc for the A3 simulation). 

2.2. Column Densities & Quasar Attenuation 

We calculate the column density between a black hole and a 
hypothetical observer from the simulation outputs spaced ev- 
ery 10 Myr before and after the merger and every 5 Myr dur- 
ing the merge r of each galaxy pair. The calculation method 
is described in Hopkins et al. (2005a), but we review it briefly 
here. We generate ~ 1000 radial lines-of-sight (rays), each 
with its origin at the black hole particle location and with di- 
rections uniformly spaced in solid angle d cos 8 d<j). For each 
ray, we then begin at the origin, calculate and record the lo- 
cal gas properties using GADGET, and then move a distance 
along the ray Ar = rjh sm i, where rj < 1 and h sm \ is the local 
SPH smoothing length. The process is repeated until a ray is 
sufficiently far from its origin (> 100 kpc). The gas properties 
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along a given ray can then be integrated to give the line-of- 
sight column density and mean metallicity. We test different 
values of 77 and find that gas properties along a ray converge 
rapidly and change smoothly for r\ = 0.5 and smaller. We sim- 
ilarly test different numbers of rays and find that the distribu- 
tion of line-of-sight properties converges for > 100 rays. 

Given the local gas properties, we use the GADGET multi- 
phase model of the ISM described in Snri ngel & Her nauist 
(2003) to calculate the local mass fraction in "hot" (dif- 
fuse) and "cold" (molecular and HI cloud core) phases of 
dense gas and, assuming pressure equilibrium between the 
two phases, we obtain the local density of the hot and cold 
phase gas and the corresponding volume filling factors. The 
values obtained correspo nd roughly to the fiducial values of 
iMcKee & Ostrikerl ill 9771) . Given a temperature for the warm, 
partially ionized component ~ 8000 K, determined by pres- 
sure equilibrium, we further calculate the neutral fraction of 
this gas, typically ~ 0.3-0.5. We denote the neutral and total 
column densities as Nm and A^h, respectively. Using only the 
hot-phase density allows us to place an effective lower limit 
on the column density along a particular line of sight, as it as- 
sumes a ray passes only through the diffuse ISM, with > 90% 
of the mass of the dense ISM concentrated in cold-phase 
"clumps." Given the small volume filling factor (< 0.01) and 
cross section of such clouds, we expect that the majority of 
sightlines will pass only through the "hot-phase" component. 

Using Lboi = £,Mc 2 , we model the for m of the in- 
trinsic quasar continuum SED following iMarconi et al.l 
J2004V base d on optical t hrough h ard X-ray obs ervations 
(e.g..lElvis et alJl994llGeorge et alJl998UYanden Berk et al] 
| 200lUPerola et alJl2002b iTelfer et al. 112002b lUeda et alJl2003t 
Vigna li et alJ 1200311 This gives a B-band luminosity 
log(Ls) = 0.80-0.067£ + 0.017£ 2 -0.0023£ 3 , where £ = 
log(Lboi/i©) - 12, and we take \g = 4400 A. We then use a 
gas-to-dust ratio to determine the extinction along a given line 
of sight at this frequency. Observations suggest that the ma- 
jority of the population of reddened quasars have reddening 
curves similar to that of the Small Magellenic Cloud (SMC; 
Hopkins et al. 2004), which has a gas-to-dust ratio lower than 
the Milky Way by approximately the same factor as its metal- 
licity (Bouchet et al. 1985). We therefore consider both a gas- 
to-dust ratio equal to that of the Milky Way, (A b /Nhi)mw = 
8.47 x 10~ 22 cm 2 , and a gas-to-dust ratio scaled by metallic- 
ity, As/Nm = (Z/0.02)(Afl/jVHT )Mw- For b oth cases we use 
the SMC-like reddening curve of Pea (I 19921) . We calculate ex- 
tinction in X-ray frequencies (0.03- 10 keV) using the photo- 
electri c absorption cross sections of Morrison & McCammon 
i ll 9831) and non-relativistic Compton scattering cross sections, 
similarly scaled by metallicity. In determining the column 
density for photoelectric X-ray absorption, we ignore the in- 
ferred ionized fraction of the gas, as it is expected that the 
inner-shell electrons which dominate the photoelectric ab- 
sorption edges will be unaffected in the temperature ranges 
of interest. We do not perform a full radiative transfer calcu- 
lation, and therefore do not model scattering or re-processing 
of radiation by dust in the infrared. 

2.3. Quasar Lifetimes & Sensitivity to Simulation Parameters 

We q uantify the quasar lifetime fg as in iHopkins et al.l 
(2005a), as a function of the limiting B-band luminosity 
^B,min- For each sightline, an observed lifetime is determined 
as the integrated time in the simulation during which the given 
sightline sees a B-band luminosity above a given threshold, 
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FIG. 1. — Quasar lifetimes as a function of L B mm from the A3 simula- 
tion. The dotted lines show the lifetime if attenuation is ignored (time with 
observed B-band luminosity Ls^bs > Lb .min)- The thick solid lines indicate 
the median observed lifetime, with thin lines giving the 25%-75% contours. 
Lifetimes are calculated using: (A) cold-phase ISM density, (B) hot-phase 
density (ignoring metallicity and ionization), (C) hot-phase density (includ- 
ing metallicity and ionization), (D) hot-phase density with artificially low 
metallicity (0.1 times solar). 



^fi.obs > ^B.min- We consider both the intrinsic lifetime tg, the 
total time the intrinsic L B > L B ^ m (ignoring attenuation), and 
the observed lifetime derived from a particular column density 
calculation. We note that below 5 Myr (the simulation output 
frequency in the standard case) our estimates of fg become 
uncertain owing to the effects of quasar variability and our 
inability to resolve the local small-scale physics of the ISM. 
However, the majority of sightlines see lifetimes in the range 
10 -20 Myr up to L^min = Lq, in good agreement with 
observations suggesting lifetimes ~ 10 7 yr and well above this 
limit. 

Figure ^ shows the lifetime obtained in the A3 simulation 
using different calculations of the column density Nm (cases 
A-D). The uppermost curve in all cases is the intrinsic life- 
time t l Q. The lower thick curve is the median observed lifetime 
calculated using a particular estimate of the column density, 
with the thin curves representing the 25%-75% inclusion con- 
tours. The cold-phase ISM density (A) is large enough for the 
quasar to be completely extincted out of observable ranges for 
the duration of the quasar phase, rendering the object observ- 
able only during quiescent phases with L B ~ 10 9 L Q . We also 
show the significantly longer lifetimes calculated using the 
hot-phase ISM density, first ignoring corrections for metal- 
licity and ionization (B; i.e. assuming A B /Nm = (A b /Nhi)mw 
and Nm = A^h) and then including these effects (C), which in- 
creases the median lifetime by a factor ~ 2 and significantly 
reduces the scatter towards shorter lifetimes. Finally, we cal- 
culate the lifetime using the hot-phase density and an SMC- 
like gas-to-dust ratio, (A B /N m )sMC = 0.869 x 10" 22 cm 2 (D; 
essentially assuming a metallicity 0.1 times solar), which sets 
a strong upper limit on the observed lifetime. These lifetimes 
are still only ~ 2 times as long as the lifetimes calculated 
using metallicity-weighted neutral column densities, and still 
~ 1/3 the intrinsic lifetime at Ls.mm = 10 11 L©. We therefore 
expect that, after accounting for the clumping of most mass 
in the cold phase of the most dense regions of the ISM, the 
qualitative relation between our observed and intrinsic calcu- 
lated lifetimes should be relatively insensitive to the details 
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FIG. 2. — Quasar lifetimes as a function of Lg m j„ from the A3 simula- 
tion, showing the waveband dependence of the lifetime. The thick lines give 
the lifetime if attenuation is ignored, thin lines indicate the median observed 
lifetime (calculated using the metallicity-weighted, hot-phase neutral column 
density). Upper panel: From black to blue to red, lifetimes are calculated at 
3000, 4500, 6000, 7500, and 9000 A. Lower panel: lifetimes are calculated 
for soft (0.5 -2keV, red) and hard (2- lOkeV, black) X-ray bands. 



of the column density calculation, with variation by a factor 
of ^2-3 between lifetimes calculated using different pre- 
scriptions for the multiphase ISM, ionization, and metallicity 
effects. 

3. FREQUENCY DEPENDENCE OF THE QUASAR LIFETIME 

We expect the intrinsic quasar lifetime above a given lumi- 
nosity in a particular band to change with the frequency of that 
band. The shift in lifetime with frequency is approximately 
given by an offset in luminosity corresponding to the dif- 
ference in luminosities at different frequencies in our model 
quasar spectrum. However, this change is not entirely history- 
independent, as the model spectrum alters shape with varying 
luminosity owing to bolometric corrections. More important, 
our model predicts that the observed lifetime above a given 
luminosity will change with frequency not only as a result of 
the shift in intrinsic lifetime, but primarily as a result of vary- 
ing levels of attenuation at different frequencies. 

Figure |2] shows the lifetimes as a function of luminosity 
for the A3 simulation, for various frequencies in the opti- 
cal and X-ray bands. We plot lifetimes in the manner of 
Figure [l] i.e. the total time with vL v > L vm i„ (L u ,min = 
(vdL/dv) m i„), at the representative optical wavelengths A = 
3000, 4500, 6000, 7500, and 9000 A (colored black to blue to 
red, respectively). We find that the intrinsic lifetime systemat- 
ically decreases at longer wavelengths (based on the shape of 
the quasar SED), while the observed lifetime increases (ow- 
ing to weakening attenuation), with the two approaching one 
another near infrared wavelengths. At longer wavelengths, we 
expect the curves to cross, as the observed lifetime in the far 
IR, owing to dust re-processing of quasar radiation, will be 
longer than the intrinsic lifetime, but we do not yet incorpo- 
rate such re-processing into our column density calculations. 
At very low L v .„„■„, lifetimes approximately converge to the 
total duration of the simulations, as the quasar is above such 
luminosities throughout almost the entire merger. 

The lower panel of Figure |2] shows the lifetimes in the soft 
and hard X-ray bands, defined as 0.5-2keV and 2- lOkeV, 
respectively. We find that the soft X-ray band is heavily at- 
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FIG. 3. — Distribution of "hot-phase", neutral column densities Nm from 
the A3 simulation. The probability of seeing a given Nu i is weighted by the 
time such an Nm is observed along all sightlines. Plots are scaled linearly 
(left) and logarithmically (right). The top panel shows all Nm values from 
all times, the lower panels give the distribution for times where the observed 
B-band luminosity is above the Lflnun shown. The smooth curve is the distri- 
bution of Nm values from the observed SDSS quasar sample ofHo pkins et alj 
i 2004), who used an absolute ;'-band magnitude limit approximately equiva- 
lent to Z-B.obs > 10 11 Z.q, scaled to the appropriate modal Nm- 



tenuated and has much shorter observed than intrinsic life- 
times, primarily owing to photoelectric absorption. However, 
the hard X-ray band is relatively unaffected, as the hot-phase 
column densities tend to be well below the Compton-thick 
regime A^h ~ 10 24 cm" 2 . Thus, most quasars obscured in the 
optical for much of their high accretion rate lifetimes should 
be observable in hard X-rays. We defer a full treatment of the 
difference between the attenuation of soft and hard X-rays to a 
future paper, but note that this difference may account for the 
slope of the cosmic X-ray background, with a population of 
obscured quasars at spectral energies < a few keV as a natural 
stage in the evolution of quasar activity. 

4. TYPICAL COLUMN DENSITIES 

Given a calculation of the column density along multiple 
lines-of-sight to the simulated quasars, we can examine the 
typical column density distributions at all times and at times in 
which the quasar is observed above some luminosity thresh- 
old. Following the analysis in Hopkins et al. (2005a), we ex- 
pect the characteristic column densities to rise rapidly during 
the merger, as gas inflow traces and fuels a rising accretion 
rate. Eventually, column densities rapidly fall in a blowout 
phase once the black hole has reached a critical mass, creating 
a window with rapidly changing column densities in which 
the quasar is observable until accretion rates drop below those 
necessary to fuel qu asar activity. As is clear in Figure 2 of 
Hopki ns et alJ fc005al) . the dispersion in A^h at any particular 
time in the simulation is generally small, a factor of ~ 2 in ei- 
ther direction, but typical A^h values can change by an order of 
magnitude over timescales ~ Myr. Therefore, it is of interest 
to consider the probability of an observer, viewing the object 
at a random time, seeing a given column density. 

Figure[3]plots the distribution of column densities Nm cal- 
culated using the hot-phase metallicity-weighted ISM density 
in the A3 simulation. The probability of seeing a given Nm 
is calculated as proportional to the total time along all sight- 
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lines that such an Nm is observed. The plots are scaled both 
linearly (left panels) and logarithmically (right panels) in col- 
umn density. Typical column densities are distributed about 
2Vhi = 10 21 cm -2 , approximately symmetrically in logiVei- We 
further plot the distribution of Nm values requiring that the 
observed B-band luminosity be above some threshold Ls.mm- 
The curve shown in the linear plots is the best-fit to the 
En-v distribution of b right SDSS quasars with z < 2.2, from 
Hopkins et al. (2004). The curve has been rescaled in terms 
of the column density (inverting our gas-to-dust prescription) 
and plotted about a peak (mode) A^hi (undetermined in Hop- 
kins et al. 2004) of N m ~ 0.5 x 10 2I cm" 2 . The /-band 
absolute magnitude limit imposed in the observed sample, 
M, < -22, corresponds approximately to our plotted B-band 
limit L^ obs > 10 n L Q . 

The agreement between the observed column density dis- 
tribution and the result of our simulations once the same se- 
lection effect is applied is strong evidence in favor of our 
model for quasar evolution. The distribution as a function 
of limiting luminosity is a natural consequence of the dynam- 
ics of the quasar activity. Throughout much of the duration 
of bright quasar activity, column densities rise to high lev- 
els as a result of the same process that feeds accretion, natu- 
rally producing a reddened population of quasars (a red "tail" 
in the quasar color distribution), extending to very bright 
quasars strongly reddened by large Nm- The existence of 
this extended reddened quasar populati on in radio and opti- 
cally selected quasars is well known (e . g.lW ebster et aljl995t 
Brothe^ri^y^y_|20 l 0Jl iFigncj^yriJ |20 01l iRichardsetall 
l200lHGregg et al.ll2002HWhite et alJl2003l) . and we are able 
to reproduce its distribution not as a distribution of source 
properties but as a result of the time evolution of quasar phe- 
nomena. Our estimate of the distribution of A^hi for Lg.obs > 
lO n L0 does not account for an additional selection effect, 
namely that strongly reddened quasars may not be extincted 
from an observed sample (if the intrinsic luminosity is large 
enough), but their colors may be significantly reddened to the 
point where color selection criteria of quasar surveys will not 
include them. However, this effect would only serve to bring 
our distribution into better agreement with observations, as it 
would slightly lower the high-A^i tail. 

Moreover, our estimates of the distribution of Nm allow us 
to make a directly observable prediction, namely that the col- 
umn density distribution in optically selected samples with a 
minimum absolute magnitude or luminosity should broaden 
in both directions (to larger and smaller Nm) as the limiting 
selection luminosity is decreased (limiting absolute magni- 
tude increased). This is because, at lower luminosities, ob- 
servers will see both intrinsically bright periods extincted by 
larger column densities (broadening the distribution to larger 
Nm values) and intrinsically faint periods with small col- 
umn densities (broadening the distribution to smaller Nm val- 
ues). The distribution of A^h values in hard X-ray quasar sam- 
ples, which are much less affected by this extinction (see §3), 
should also reflect this trend, and indeed the distribution of 
A^h values is much broader and flatter, with quasars seen from 
Ah ~ 10 20 - 10 24 crrT 2 dUeda et al.l2003l) . 

Extension of this distribution to the other black hole masses 
and peak quasar luminosities of our simulations reveals sim- 
ilar qualitative behavior. The characteristic column densities 
during the obscured, intrinsically bright phases of accretion 
generally increase with galaxy and black hole mass, but typ- 
ically by factors < 10 across the range of our simulations 
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FIG. 4. — Upper panel: Thick contours plot the median column density 
as a function of simulation time, with thin contours at 25% and 75% inclusion 
levels, for the A4 simulation (V v j r = 226 kms~'). Solid contours represent the 
density of the "hot-phase" ISM, dashed contours the total simulation density. 
Middle panel: Bolometric luminosity of the black hole in the simulation, 
Ami = trMc 1 (thick), and ratio of bolometric to Eddington luminosity, / = 
Lh a \/LEdd (thin). Values are shown for each simulation timestep. Lower 
panel: Observed B-band luminosity with Nm calculated using the median 
"hot-phase" ISM density. 



V v ir = 80 -320 km s . Additionally, although mean column 
densities increase somewhat with increasing mass, the func- 
tional form of the Nm distribution remains similar at each 
limiting luminosity, implying that a combined population of 
quasars similar to those we have simulated will match the ob- 
served Nm distribution with appropriate selection effects ac- 
counted for as in Figure [3] Finally, we note that although we 
do not see extremely large column densities A^h > 10 cm" 2 
in the distributions from our simulations, our model does not 
rule out such values. It is possible that very bright quasars 
in unusually massive galaxies or quasars in higher-redshift, 
compact galaxies which we have not simulated may, during 
peak accretion periods, reach such values in their typical col- 
umn densities. Moreover, as our model assumes ~ 90% of the 
mass of the densest gas is clumped into cold-phase molecular 
clouds, a small fraction of sightlines will pass through such 
clouds and encounter column densities similar to those shown 
for the cold phase in Figure^ as large as A^h ~ 10 26 cm" 2 . This 
also allows a large concentration of mass in sub-resolution 
obscuring structures, such as an obscuring toroid on scales 
< lOOpc, although many of the phenomena such structures 
are invoked to explain can be accounted for through our model 
of time-dependent obscuration. 

5. EXTENSION OF QUASAR LIFETIMES TO OTHER MASSES 

We generalize the results from iHopkins et alJ 1120051 
and those above (which use simulation A3) to a range of 
masses, using the simulations Al, A2, A4, and A5. Re- 
call, these simulations are identical except for increasing 
galaxy mass given by the sequence in virial velocity, V vlr = 
80, 113, 160,226, and320kms" 1 for A1-A5, respectively. The 
evolution of these simulations would be simply related if not 
for the scale-dependent physics of cooling, star formation, 
and black hole accretion and feedback. Figures @] and [5] plot 
the bolometric luminosity of and column densities to the su- 
permassive black hole as a function of simulation time, for the 
period shortly before and after the merger in simulations A4 
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FIG. 5. — Same as Figure|4] but for the A5 (V v ii- = 320 kms ') simulation. 



and A5, respectively. This is presented in the style of Figure 
2 from Hopkins et al. (2005a), which showed the same for 
A3. The median column densities and 25%-75% contours are 
indicated for the total average simulation density, and the cal- 
culated hot-phase density, as described in §2.2, along with the 
observed B-band luminosity using the median hot-phase A^hi- 
The final masses of the black holes are also shown, where 
M BH ,finai ~ 7 x 10 6 ,3 x 10 7 ,3 x 10 8 ,7 x 10 8 ,and 2 x 10 9 for 
simulations A1-A5, respectively. The evolution seen in sim- 
ulation A4 is quite similar to that of A3, with a sharp rise at 
the time of the merger to very large peak Eddington-limited 
luminosities, for a period ~ 10 8 yr, much of which is unob- 
servable in the optical owing to a corresponding rise in Nm- 
At the final stages of strong accretion, when gas is rapidly 
being consumed and expelled, the observable luminosity fi- 
nally rises to peak levels over a much shorter time interval 
~ 10 7 yr. In the A5 simulation, the period of rapid accretion 
and large, Eddington-limited intrinsic luminosities is signifi- 
cantly longer, ~3x 10 8 yr, and the peak luminosity is higher, 
as expected. However, as before, most of this period is ob- 
scured by large column densities, and, remarkably, once again 
only the final stages of intrinsically high luminosity are visi- 
ble, as gas is expelled and the accretion rate and column den- 
sity both fall. The same behavior is seen in the Al and A2 
cases, though with much lower peak bolometric luminosities 
~ 10 10 and 1O H L0, respectively. The striking similarity of 
this process across all our simulations suggests that it is at 
least qualitatively insensitive to the details of the accretion 
history and galaxy or black hole masses. 

We quantify the resulting quasar lifetimes in Figure [6j 
where we plot the quasar lifetime tQ as a function of the lim- 
iting B-band luminosity LB, mm - The panels show the lifetimes 
for simulations A1-A5, from top to bottom. In each panel, the 
solid curve shows the intrinsic lifetime t Q ; i.e. the total time 
the intrinsic Lb > ifl.min- The dashed curves show the inte- 
grated time that the observed B-band luminosity meets this 
criterion, using the median metallicity-weighted hot-phase 
neutral column density of the simulat ion (case C of F igurelD. 
As demonstrated in Figure ^ and in Hopki nset alJ (|2005a), 
varying this definition of the column density over realistic 
ranges can change the intrinsic lifetime by a factor of ~ 2-3 at 
a given LB iimn ■ We also plot the lifetime as a function of Z, the 
ratio of the bolometric luminosity to the Eddington luminos- 
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FIG. 6. — Left panels: Quasar lifetimes as a function of Ls, m i n from the 
simulations (time with observed B-band luminosity L s i, s > L B m i n ). Results 
from simulations A1-A5 are shown from top to bottom. Solid lines indicate 
the lifetime if attenuation is ignored. The dashed lines give the median ob- 
served lifetime calculated with the metallicity-weighted, hot-phase neutral 
column density. Right panels: Same as the left panels, but with fg defined 
as a bolometric luminosity Lb i > IL^j, where the observable lifetime is the 
time above a given ratio I with B-band optical depth less than unity. 

ity, and find a similar trend. Here, we calculate the "observed" 
time above a given ratio / as the time above such a ratio with 
a B-band optical depth less than unity. In all cases, the ob- 
served tQ is significantly smaller than the intrinsic lifetime for 
all Lb. min > 10 9 Lq, and the ratio fg/fg decreases with increas- 
ing Lfi.min. However, both the intrinsic and observed lifetimes 
increase systematically at a given Lb > 10 9 L Q with increas- 
ing galaxy and black hole mass, and the peak luminosity sim- 
ilarly increases, as expected. In all cases corresponding to 
quasar-like luminosities, our estimates of observed lifetimes 
of ~ 10 7 yr agree well with observations, with much longer 
intrinsic lifetimes > 10 8 yr. 

A quasar accreting and radiating at the Eddington limit, 
with luminosity increasing exponentially up to some peak lu- 
minosity L max after which accretion shuts off, will spend an 
integrated time above a given luminosity L given by 

t Q (L' >L) = t s \n(L max /L), (1) 

where t s is the Salpeter time, t s = M B u/M = 4.5 x 10 7 yr at the 
Eddington rate with our efficiency e r = 0.1, and L max is given 
by the Eddington luminosity of the final black hole mass. This 
lifetime agrees well with the plotted lifetimes for simulations 
A3, A4, and A5 at the high-luminosity end. However, sim- 
ulations Al and A2 do not radiate at their Eddington limits 
for a long period of time, and the slope of this relation is too 
shallow in all cases, underpredicting the time that the simula- 
tion spends at luminosities L < 0.1L max . This is because the 
black holes in the simulations spend significant time at lower 
Eddington ratios both going into and coming out of the quasar 
stage and in an extended quiescent phase. Moreover, the ob- 
served lifetimes do not at all correspond to the lifetimes pre- 
dicted by this simple model. Therefore, to simply describe the 
intrinsic and observed lifetimes as a function of limiting lumi- 
nosity and final black hole mass (or peak luminosity L max ), we 
fit the calculated lifetimes tQ to truncated power laws, from 
L = 1O 9 L0 to approximately 0.5L max , or where tQ < lMyr, 
whichever occurs first. Thus we fit to the form 

tQ(L'>L) = t 9 (L/L 9 ) a , (2) 
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FIG. 7. — The results of fitting quasar lifetimes tQ to truncated power laws 
of the form 1q(L' > L) oc L a . Upper panel shows the intrinsic (thick solid 
line) and observed (thin solid line) lifetimes in the A3 simulation, with dashed 
lines indicating the corresponding best-fit truncated power laws. The dotted 
line shows the lifetime expected for Eddington-limited accretion up to the 
final black hole mass. The lower panels plot the power law slopes a and 
normalizations at L = 10 9 Z.q of the power law fits to the intrinsic (thick) and 
observed (thin) lifetimes from the series of simulations A1-A5, as a function 
of the final black hole mass in the simulation. The dotted line in the left panel 
gives a simple linear fit to a as a function of the logarithm of the final black 
hole mass. 



where Lg = 10 9 L© and tg = tQ(L' > Lg). 

Figure [7] shows the resulting power-law slopes a and nor- 
malizations tg of the lifetimes fg as a function of the final 
black hole mass (or equivalently, the peak quasar luminos- 
ity given as the Eddington luminosity of the final black hole 
mass). The calculated, fitted, and Eddington-limit approx- 
imated lifetimes are also shown for comparison for the A3 
simulation. The lifetimes are well-fitted by power laws in 
this range, and the discrepancy between the power-law and 
Eddington-limit approximated lifetimes at lower luminosities 
is clear. The slope of both the intrinsic and observed lifetimes 
becomes shallower at higher black hole masses, as the life- 
time increases for all luminosities and extends to higher peak 
luminosities, but the slope of the observed lifetime evolves 
much more rapidly than the slope of the intrinsic lifetime. 
Over the range of our simulations, the intrinsic lifetime slope 
can be approximated as a = -0.79 + 0.32 log(M BH /10 8 M Q ), 
or a = -0.95 + 0.321og(L max /10 12 L Q ). Although the fluctu- 
ation in the normalization tg is large for observed lifetimes, 
there does not appear to be a systematic trend with final black 
hole mass, which is expected as by luminosities of this level or 
lower, the lifetime is dominated by the quiescent-phase life- 
time throughout the duration of the simulation. This is es- 
pecially clear in the normalization of the intrinsic lifetimes, 
which also show no systematic trend and are all within a fac- 
tor of ~ two of 2 Gyr, the approximate total duration of the 
simulations. 

In order to obtain a lower limit to quasar lifetimes at low 
luminosities, we repeat this analysis, but ignore all times 
prior to the intrinsically bright quasar phase (essentially the 
final merger). This excludes the luminosity of the black hole 
as it grows from its small seed mass, and ignores the peri- 
ods of accretion with high gas fractions early in the merger, 
giving a lower limit corresponding to an already-large black 
hole suddenly "turning on" in a quasar phase. As expected, 



the Eddington-limited phase of the lifetime is more domi- 
nant in this case. However, the pure Eddington-limited model 
still gives lifetimes too short by an order of magnitude at 
low luminosities. Fitting to power laws we find very simi- 
lar relations to those shown in Figure with slightly shal- 
lower slopes, a = -0.61 + 0. 30 log(M BH /lO 8 M ) for the in- 
trinsic lifetime, and normalization tg ~ 1 Gyr. We further fit 
one additional functional form to the intrinsic lifetimes, in 
which the lifetime follows the Eddington-limited shape for 
large luminosities, down to some luminosity Lb rea it, and be- 
low this luminosity follows a power-law. This allows a con- 
tinuous description of the lifetime, very accurate at both low 
and near-peak luminosities. For these fits, we find similar 
power-law slopes a = -0.81 + 0.32 log(M B H/10 8 M Q ), with 
Lbreak ~ 0.025 L\, i, max ~ 0.14L B ,„ MA . nearly constant. We also 
fit the lifetimes as a function of frequency for the represen- 
tative visual wavelengths shown in Figure [2] Although, as 
expected from the figure, the observed lifetime slopes be- 
come slightly shallower and the intrinsic lifetime slopes be- 
come slightly steeper as the quasar is observed at redder wave- 
lengths, the change with wavelength across this range is very 
small compared to the change with black hole mass. More- 
over, the change in normalization tg across visual wavelengths 
is consistent with zero. We therefore find that quasar lifetimes 
can be well approximated by simple power laws, where the 
power law slope depends on the final black hole mass or peak 
luminosity of a system, and find significantly longer lifetimes 
at low luminosities than predicted by simple models of ac- 
cretion at a constant Eddington ratio. Furthermore, attempts 
to de-convolve intrinsic quasar properties and observations as 
well as semi-analytical attempts to model observable quasar 
luminosity functions and statistics require an accurate model 
of quasar lifetimes as a function of luminosity and host sys- 
tem properties. The fits we describe have the advantage of a 
simple, analytical form which can be applied to a continuum 
of luminosities and intrinsic quasar properties, while address- 
ing the shortcomings of simpler lifetime models below peak 
system luminosities. 

6. DISCUSSION & CONCLUSIONS 

In this p aper, we have studied implications of the con- 
clusion of Hop kins et alJ J2005al) that incorporating the ef- 
fects of obscuration in a galaxy merger simulation gives ob- 
served quasar lifetimes of ~ 10 7 yr. We extended our analy- 
sis to a series of five simulations with virial velocities V wn - = 
80- 320 km s" 1 , producing a range of final black hole masses 
from ~7xl0 6 -2x 1O 9 M , and find qualitatively similar 
results in all cases and for different determinations of the 
obscuring column density. The good agreement with obser- 
vations and the significant difference between observed and 
much longer intrinsic lifetimes is found in all our simulations. 
The qualitative evolution of quasars in this model is robust, 
with the result reproduced across a range of galaxy masses, 
and using different methods for calculating column densities 
that can give factors greater than an order of magnitude dif- 
ferent attenuations along any given line of sight at a particular 
time. We therefore expect that quasars will have extended in- 
trinsic lifetimes, much of which are obscured in visual wave- 
lengths. The processes which fuel periods of quasar activity 
by channeling significant quantities of gas into the central re- 
gions of a galaxy will produce this extended intrinsic lifetime 
of > 10 8 yr. However, for a significant fraction of the intrin- 
sic lifetime for quasar activity, this same process will produce 
large column densities which heavily obscure the quasar and 
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attenuate it well below observable limits in the B-band and 
other visual wavelengths. Eventually, feedback from accre- 
tion energy will remove surrounding gas, creating a window 
in which the quasar is optically observable before accretion 
rates drop below those needed to maintain quasar luminosi- 
ties. 

Of course, not all aspects of AGN activity are related to 
mergers. For example, some low redshift quasars (e.g. Bah- 
call et al. 1996) and many Seyferts appear to reside in rel- 
atively ordinary galaxies. Gas expelled by quasar feedback 
may cool and relax to the central regions of the galaxy, and 
passive stellar evolution can produce a large fuel supply in 
winds. In our picture, it is most natural to interpret occur- 
rences of AGN activity that are not caused by mergers as 
arising from normal galaxies being re-activated as gas is spo- 
radically accreted by a black hole created earlier in a bright 
quasar phase. The principle requirement of this modeling is 
that such activity should not contribute a large fraction of the 
black hole mass, to avoid spoiling tight correlations between 
the black hole and host galaxy properties or overproducing 
the present-day density of very massive black holes. Ra- 
diative heating from the final stages of quasar activity may 
suppress much of this activity (much of the remaining gas 
is heated to very high temperature s; e.g. Cox et al. 2005), 
and the hydrodynami cal modeling of Ciotti & Ostriker; ( 119971 
2001): lSazonov etaflll200 5) demonstrates that residual accre- 
tion from stellar feedback and "cooling flows" is dramatically 
suppressed by even relatively small AGN feedback once the 
black hole has reached a critical mass determined by the Mbh- 
<j relation. 

Furt her, demographic arguments such as that of fSoltan 
as well as more detailed analysis by e.g., 
lYu & Tremaind ll2002l) implies that most of the present black 
hole mass density is accumulated in bright quasar phases, for 
which our modeling should be most applicable. This remains 
true in our analysis, and we demonstrate in Hopkins et al. 
(2005b,c) that our modeling can account for the entire quasar 
luminosity function in different frequencies, although we de- 
fer a detailed prediction of the resulting supermassive black 
hole density and mass function to a future paper. It is impor- 
tant to note, however, that given our prediction of an obscured 
accretion phase with a duration up to ~ 10 times that of the 
optica lly observable phase, a naive application of the Soltan 
( 1982) or other demographic arguments would imply a simi- 
lar increase in the black hole mass density accumulated in ob- 
scured quasar phases. But this neglects both the fact that black 
holes are growing during this phase and the luminosity de- 
pendence of the quasar lifetime. A more detailed calculation 
shows that a large fraction (> 80% in the brightest quasars) 
of the quasar mass is accumulated in the optically observable 
stages of black hole growth when the final e-folding of expo- 
nential black hole growth occurs as the black hole begins to 
drive strong feedback and soon shuts down accretion. Simi- 
larly, most of the total black hole radiated energy is emitted 
in this stage, with a considerable fraction (~ 0.3 - 0.6) ob- 
servable. Thus, we expect that corrections to the estimates of 
black hole growth from demographic arguments based on ob- 
servations of bright, relatively unobscured quasars will in fact 
only be of order unity from our modeling, and thus should not 
conflict with theoretical limits. It is still an interesting and 
important question, however, to determine from our model- 
ing the implied relic AGN density and constraints from the 
combination of observed luminosity functions and black hole 
mass functions. 



We do not expect other considerations, such as the orbital 
parameters of the merger or the collimation of black hole 
feedback, to substantially change our results. Different or- 
bital separations, orientations, and energies will change the 
time of interaction, but ultimately our model for quasar evo- 
lution depends only on the merger eventually depositing suf- 
ficient gas into the central regions of the remnant to fuel rapid 
black hole growth. After such a process begins, the evolution 
of the black hole should be determined by the self -regulated 
m echanisms which termi nate black hole growth, as described 
in lDi M atteo et al. (2005). This will eventually expel gas, ren- 
dering the quasar observable and shutting down the accretion 
phase at the final black hole mass given by the Mbh-c rela- 
tion. The phase of rapid black hole growth and the expulsion 
of gas in the final stages of strong accretion occur during the 
end-stages of the merger and in the very center of the merg- 
ing cores, when and where the structure of the original orbits 
should be least important. This view is supported by addi- 
ti onal s imulations not pre sented here but shown in e.g. Fig. 16 
of Springel et al. (2005b), demonstrating that the final black 
hole mass is relatively insensitive to orbital parameters. 

The exact mechanism of black hole feedback should sim- 
ilarly not dramatically change this picture, so long as a pe- 
riod of rapid black hole growth occurs, requiring large den- 
sities (which naturally generate obscuring column densities), 
and ultimately terminates via some self-regulated mechanism 
which expels gas from central regions, rendering the quasar 
optically observable in the final stages of its life. In particular, 
collimation of the black hole feedback may significantly in- 
crease the energy or momentum input on-axis, "blowing out" 
material along this axis earlier in the quasar lifetime. How- 
ever, since the black hole is growing exponentially during 
peak accretion, this should occur at most on the order of a 
couple of Salpeter times (during which the black hole lumi- 
nosity grows by an order of magnitude) before it would in 
the isotropic case, changing the intrinsic lifetime by only a 
fraction of its total. The unobscured stage of the quasar life- 
time may then begin slightly earlier along on-axis sightlines, 
but should be of comparable duration as a similar mechanism 
removes the gas from central regions. In addition, in order 
for significant black hole growth to cease at the same critical 
mass as in the isotropic case (determined by the Mbh-c re- 
lation), the high-accretion rate phase must end over the same 
timescale ~ t$ as the "lead" with which collimated, on-axis 
blowout precedes isotropic blowout. The necessity that self- 
regulation eventually terminates accretion implies that sur- 
rounding densities, even off-axis, must decrease, resulting in 
at most a broadening of the column density distribution as typ- 
ical column densities off-axis could be somewhat larger than 
column densities on-axis. 

For much of this regime, we find that the source should 
be observable in hard X-ray frequencies, as the attenuation is 
weaker and column densities generally lie below the Comp- 
ton thick limit. This distinction naturally produces a substan- 
tial population of obscured quasars, as a standard phase in the 
evolution of quasars over their lifetimes. Although we de- 
fer a detailed calculation of the quasar luminosity function as 
a function of redshift and in different observed frequencies 
(Hopkins et al. 2005b), we note that such a population can 
naturally account for both the cosmic X-ray background spec- 
trum and discrepancies between optical and X-ray luminosity 
functions at various redshifts (Hopkins et al. 2005c). More- 
over, optical and X-ray observational samples will be affected 
rather differently by selection effects and magnitude limits, 
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especially for quasars near their peak luminosities. Therefore, 
the differences in the evolution of optical and X-ray selected 
samples may be accounted for through the dependence of ob- 
scuration on intrinsic luminosity and host galaxy properties. 
Attenuation is also weaker in the infrared, and the large frac- 
tion of optical and UV energy absorbed during the obscured 
phase may be reprocessed by dust, appearing as thermal ra- 
diation in the IR. This may render the observable lifetime 
larger than the intrinsic lifetime at long wavelengths, an ef- 
fect important for calculating the IR background and estimat- 
ing source populations. Similarly, the total obscuration and 
obscuration as a function of intrinsic luminosity are impor- 
tant for the evolution of luminous and ultraluminous infrared 
galaxies, and for estimating the relative energy contributions 
in the IR of starbursts and active galactic nuclei (AGN). 

The typical column densities in our simulations corre- 
spond well to observed column densities of optically selected 
quasars, once the appropriate observed magnitude limit has 
been imposed. Thus, our model allows us to reproduce the 
distribution of optically obscured sources above this magni- 
tude limit naturally, as an evolutionary effect of the mecha- 
nisms which fuel and regulate quasar growth, without invok- 
ing particular distributions of source properties or geometric 
patterns of obscuration. Moreover, the column densities we 
calculate in our simulations allow us to predict that the distri- 
bution of jVh values should become broader as the minimum 
observable luminosity is decreased, as both faint, quiescent 
phases with low (< 10 20 cm" 2 ) column densities and bright, 
obscured phases with high (> 10 22 cm" 2 ) column densities be- 
come observable. This effect is seen in the broad distribution 
of Nh values in X-ray samples, which are much less affected 
by attenuation, but a more detailed analysis including model- 
ing the distribution of quasar properties is needed to reproduce 
both the optical and X-ray A^h distributions more accurately. 

We find that the peak luminosities and lifetime above any 
given luminosity increase systematically with galaxy mass 
and final black hole mass, although the distinction between 
observed and intrinsic lifetimes remains significant in all 
cases. Moreover, the ratio of observed to intrinsic lifetimes 
decreases in all cases with increasing minimum luminosity. 
We also find that intrinsic lifetimes at luminosities < 0. 1 times 
the peak luminosity are poorly fit by assuming quasars always 
accrete near the Eddington limit, but rather that lifetimes are 
well-fitted by power laws with a steeper slope. This is a re- 
sult of sub-Eddington accretion rates both before and after 
the peak accretion phase, even at luminosities significantly 
above the quiescent steady-state small accretion rates seen at 
late times in the simulations. This suggests that many quasars 
seen at low luminosities may be quasars with a large peak lu- 
minosity in a significantly sub-Eddington phase, although the 
fraction of quasars observed at high Eddington ratios becomes 
large with increasing luminosity. 

Additionally, these lifetimes imply that, even at high lu- 
minosities where growth may be Eddington-limited, quasars 
spend a significant fraction of their lives with intrinsic lumi- 
nosities well below their peak luminosities. Therefore, any 
observed luminosity function is the convolution of the distri- 
bution of quasars with a given peak luminosity (determined by 
the final black hole mass and thus the merging galaxy proper- 
ties) and a non-trivial light curve (Hopkins et al. 2005b). It 
is clear from these calculations that any attempt to theoreti- 
cally model even the intrinsic luminosity functions of quasars 
must take into account the functional dependence of the light 
curve on luminosity and time. In order to apply these models 



to observed luminosity functions, the dependence of observed 
luminosity on intrinsic luminosity and the quasar evolution, as 
well as the observed frequency, must be considered (Hopkins 
et al. 2005b,c). 

We find that the ratio of observable to intrinsic lifetimes is 
a strongly decreasing function of the limiting luminosity of 
observations. Modeling this effect is necessary to estimate in- 
trinsic quasar lifetimes from observations, as well as for using 
theoretically motivated accretion models to predict the quasar 
luminosity function and space density of present-day super- 
massive black holes. Further, the effects we describe can ac- 
count for the presence of an obscured population of quasars 
which are missed by optical, UV, or soft X-ray surveys but 
may contribute significantly to the cosmic X-ray background. 
Observations of the cosmic X-ray background and compari- 
son of the optical and hard X-ray quasar luminosity functions 
that indicat e the existence of a lar ge obscured population of 
quasars (Brandt & Hasinger 2005, and references therein), is 
explained in our picture, because it predicts different observed 
lifetimes and populations at different frequencies. The sce- 
nario we describe also implies that the reprocessing of quasar 
radiation by dust in surrounding regions can account for ob- 
servations of luminous and ultraluminous infrared galaxies 
with merger activity and obscured AGN. 

Together with the modeling presented by Di Matteo et 
al. (2005), Springel et al. (2005a,b), Springel & Hern- 
quist (2005), Hopkins et al. (2005a,b,c), and Robertson et 
al. (2005), the results described here motivate the follow- 
ing picture for galaxy formation and evolution. Through the 
hierarchical growth of structure in a cold dark matter uni- 
verse (White & Rees 1978), mergers between galaxies oc- 
cur on a regular basis. Those involving gas-rich progen- 
itors, which would be increasingly more common towards 
higher redshifts , produce inflows of gas through gravitational 
torques ([B arnes & Hernaui slll991iri996h . causing starbursts 
( Mihos & Hernauist 1994, 1996) like those associated with 
luminous infrared galaxies (e.g. Sanders & Mirabel 1996). 
The high gas densities triggering these starbursts fuel rapid 
black hole growth. For most of the period over which black 
hole growth occurs, optical quasar activity would be buried, 
but X-rays from the black holes explain the p resence of non- 
therm al point sources in e.g. NGC 6240 ( Komos sa et alJ 
2003), and reprocessing of most of the black hole energy 
by surrounding gas and dust can, in principle, account for 
the spectral energy distributions of "warm" ULIRGs (Sanders 
et al. 1988c) and recent observations of a correlation be- 
twee n quasar obscuration and far-infrare d host luminosity 
(e.g-lPa ge et all2004UStevens et al.l2005l) . As the black hole 
mass and radiative output increase, a critical point is reached 
where feedback energy starts to expel the gas fueling accre- 
tion. For a relatively brief period of time, the galaxy would be 
seen as an optical quasar with a B-band luminosity and life- 
time characteristic of observed quasars. This phase of evolu- 
tion is brief (~ 10 7 yr), owing to the explosive nature of the 
final stages of black hole growth as the gas responds dramat- 
ically to the feedback energy from the exponentially evolving 
black hole. This AGN feedback terminates further black hole 
growth, leaving a remnant that resembles an ordinary galaxy 
containing a dead quasar and satisfying the Mbh-c relation. 
Modeling both the dependence of quasar lifetime on lumi- 
nosity, and the complex, time-dependent evolution of quasar 
obscuration is thus crucial to any observational or theoretical 
understanding of quasars and quasar host galaxy evolution. 
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